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Abstract 

In this paper we obtain unstable even-parity eigenmodes to the static regular 
spherically symmetric solutions of the SU(2) Yang-Mills-dilaton coupled system of 
equations in 3 + 1 Minkowski space-time. The corresponding matrix Sturm-Liouville 
problem is solved numerically by means of the continuous analogue of Newton's 
method. The method, being the powerful tool for solving both boundary-value and 
Sturm-Liouville problems, is described in details. 

1 Introduction 

Regular stationary/static solutions of nonlinear field equations play a very important role 
for understanding of non-perturbative aspects of field dynamics in various modern field 
models, including gravity. During the last decade since Bartnik and McKinnon's discovery 
PP of the regular particle-like solutions to the coupled system of the Einstein- Yang-Mills 
(EYM) equations there was a growing interest in revealing their possible physical and 
mathematical significance [2j ^ |S] • At the same time, the qualitatively similar (essentially 
non-Abelian ones) regular solutions were found in other various field models, i.e. in the 
coupled system of the Yang-Mills-dilaton (YMD) equations jB], in the coupled system of 
the Einstein- Yang-Mills-dilaton (EYMD) equations j7j- jTU], in more sophisticated stringy 
inspired field models jTTj etc (see Ref. |E] for a review). 

The main physical and mathematical properties of the mentioned solutions are deter- 
mined by the nonlinear nature of the SU(2) Yang-Mills (YM) field which is a substantial 
part of all these models. All the mentioned coupled systems of equations admit spherically 
symmetric particle-like solutions in which the Yang-Mills field configuration is bounded 
by gravity (EYM), or by the scalar dilaton field (YMD), or by both of them (EYMD) 
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in unstable equilibrium. (Note, pure YM field in the 3 + 1 Minkowski space-time is 
self- repulsive and no regular pure YM configurations exist.) In terms of the linear per- 
turbation analysis these solutions have even-parity unstable modes jlSI, [Oj, which are 
called gravitating (EYM) or dilatonic (YMD) unstable modes. 

On the other hand, in the case of a spherical symmetry the SU(2) Yang-Mills self- 
interaction potential has a double-well shaped form and its two distinct minima corre- 
spond to the topologically distinct YM vacua. The relevant YM magnetic function must 
tend to the vacuum values at the origin and at the spatial infinity for all regular solu- 
tions. It was shown some those regular solutions could play a role at ultramicroscopic 
distance analogous to that of electroweak sphalerons |2] due to existence of the odd- 
parity Yang-Mills negative modes |SI (which are called sphaleronic unstable modes) and 
the corresponding fermion zero modes, which are responsible for the anomalous fermion 
production [Hj. 

However, in this paper we concentrate on study of the even-parity unstable modes in 
the YMD system, since their significance is inspired by an interesting story, related to 
the singularity formation in nonlinear evolution equations. Indeed, it has been discovered 
in the numerical studies of the massless fields collapse jHj, there are two types of the 
collapse behavior. Type-II behavior is characterized by the mass gap absence in the black 
hole spectrum hence the black holes with an arbitrary small mass can be obtained in 
this way. Later on, the Type-I behavior was observed as well JHI- It takes place 
if the considered system of Einstein-matter equations admits regular static finite energy 
asymptotically flat solutions. In this case the smallest black hole has a finite mass which 
is equal to the mass of the lowest static solution. In the EYM system the lowest static 
solution is the = 1 Bartnick-McKinnon (BK) solution and its mass determines the 
smallest mass of the formed black hole in this type of a collapse. = 1 BK solution 
is occurred to be an intermediate attractor, which the collapsing solution should attain 
in order to turn to the Type-I black hole formation scenario ^2], ^Z|, ^Hj- Moreover, 
N = 1 BK solution is occurred to be a threshold configuration, which separates the black 
hole formation and dispersion scenarios. The revealing of the threshold nature remains 
an unresolved problem and the studies of the static configurations decay via their proper 
unstable eigenmodes seem to be a very natural step towards understanding the threshold 
nature. 

Thus, the goal of our present paper is to obtain unstable eigenmodes of the regular 
static solutions in the YMD system. It has been realized recently that the singularity 
(black hole) formation in gravity and blow-up in nonlinear wave equations share many 
common features [TH] , [201 , , [22] , since one deals with the same class of the supercritical 
evolution PDE's |2S1 • It was shown in Ref . (21] , a coupled system of Yang-Mills-dilaton 
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equations in 3 + 1 Minkowski space is also a supercritical system in PDE terminology and 
it does exhibit blow-up (singularity formation) behavior, which is, in principle, equivalent 
to the black hole formation. So, we study the unstable modes of the static solutions in the 
YMD system, which is representative for a wide class of a supercritical evolution PDE's, 
including self-gravitating systems. 

Our main tool to attack the problem numerically is a continuous analogue of the 
Newton's method (CANM), which is a really powerful tool for solving both boundary- 
value and Sturm-Liouville problems. The key point of the method is an introduction of 
formal evolution parameter t and a corresponding reformulation of the problem in terms 
of a new system of PDE evolution equations. If necessary conditions are satisfied, the 
solution of the evolution problem approaches the solution of the initial (boundary-value 
or Sturm-Liouville) ODE problem as t ^ +00. 

The paper is organized as follows. In the next Section we briefly remind the main 
features of the regular static spherically symmetric solutions to the coupled system of the 
SU(2) Yang-Mills-Dilaton equations. In the third Section the iterative scheme, based on 
the continuous analogue of Newton's method (CANM) is introduced. The CANM method 
is used for numerical solution of the matrix Sturm-Liouville problem, which describes 
(unstable) eigenmodes of the static regular YMD solutions. The fourth Section contains 
the results of the numerical experiments — the eigenvalues and the eigenfunctions of the 
mentioned Sturm-Liouville problem. We conclude with some remarks. A more detailed 
description of the CANM method is presented in the Appendix. 



2 Basic equations. 

Static spherically symmetric regular solutions of a coupled system of Yang-Mills- 
dilaton (YMd) equations in 3 + 1 Minkowki space-time were obtained and analyzed in 
Ref. jH]. In this Section we briefly remind their results. 
A coupled system of YMd fields is given by the action: 

^ = s / (5 - ^^^°'"^;^) *. (1) 

where $ is the dilaton field, F"'^'^ - the Yang-Mills field, k and g are dilaton and gauge 
coupling constants, respectively. Purely magnetic ansatz for the static SU (2) spherically 
symmetric YM potential has a form: 

k 

At = At = ea^k^[fir)-1], (2) 

hence the desired static solutions are described in terms of two independent functions: 
/(r) - the YM function and $(r) - the dilaton function. After substitution into 
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and rescaling $ — > ^/k, r — > {k/g)r, t {k/g)t, and S ^ g * kS, the dependence on 
two parameters k and g is effectively vanished. The integration over the angular variables 
gives the reduced action: 



S 



r^$'^ + exp{$} 



2r2 



dr , 



(3) 



where prime stands for derivatives with respect to the radial variable r. The corresponding 
field equations 

f{l-P) 



2$' 



2r2 



(4a) 
(4b) 



are obtained by a variation of the action (jS)) in respect to the functions /(r) and $(r) 
under the following boundary conditions: 



/(0) = ±1, /(oo) = ±l, 
$'(0)=0, <l>'(oo) = 0. 



(5a) 
(5b) 



The invariance of the equations (0]) in respect to the replacement / ^ — / allows one to 
put /(O) = 1 without loss of generality. The transformation $^$ + A, r^r exp{— A/2} 
with A = const does not change the equations either, hence one can put $o = 0. 

The system of nonlinear equations (j3)) has two singular points r = and r = oo and 
the solutions are supposed to be bounded at these points. 

Regular solutions at the point r = admit the following simple series expansions: 



/(r),_o = l-&r2 + 0(0, 
<l>(r),^o = $0 + + O(r^) 



(6a) 
(6b) 



Here 6 is a free parameter, that finally provides the boundary conditions at the origin as: 



/(0) = 1, /'(0)=0, (7a) 
$(0)=0, <l>'(0)=0. (7b) 

The asymptotical behavior of regular solutions at r — oo has the following form: 

r 
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Fig. 1: Static solutions with N = 1,2,3,4: YM function /(r) — on the left, dilaton function 
$(r) on the right. 

where c, d, and $oo are parameters. As a result, the corresponding boundary conditions 
at the infinity are 

hm /(r) = (-1)^, hm/'(r) = 0, 

r— >oo r^oo 

hm $(r) = $oo, hm <l>'(r) = 0, (9) 

r— >oo r^oo 

for Ar = 1,2, ... 

It is well-known fact that the system of equations (jH) has an infinite set of regular 
solutions, labeled by the number of the nodes of YM function /(r) [ni,|2n]- In the Ref. jB] 
the solutions have been found numerically by use of the shooting strategy in respect to the 
free parameter b. Indeed, in order to get solutions it is possible to transform the boundary 
value problem (j3)), ()7a|l . into the initial value problem with the initial conditions (jHI) 
so that the desired solutions meet the asymptotics at infinity. For the selection of the 
parameter b the following key property |3] of the YM function was used: YM function /(r) 
is strictly bounded within the interval —1 < /(r) < 1 for all globally regular solutions of 
the system (jl)). 

For our present purposes the solutions of the system (jl)) have been found using the 
continuous analogue of Newton's method jSUI (see below Chapter 3). In this approach the 
problem is considered as a boundary value problem within some finite interval r G [0, Roo], 
where < -Rqo < oo is an "actual infinity". The boundary conditions can be written as 
follows: 

/(0) = 1, /(i?oo) = (-l)^, 
$'(0) = 0, $'(i?oo) = 0. 



For numerical realization it is effective to use a boundary condition of a mixed type 
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[f{R^)-(-lf]+f'{R^) = 0, (11) 



which directly follows from the asymptotics 

The static solutions corresponding to = 1, 2, 3, 4 are shown in Fig. ^ 
Following the lines of Ref. [0] we shall consider small spherically symmetric perturba- 
tions of the obtained regular static solutions /Ar(r), $Ar(r) of the following form: 

f{r) = fN{r) + e exp{-$jv/2} v{r) exp{tujt} , (12a) 

$(r) = $7v(r) + e — u{r) exp{iujt}. (12b) 

Here e is a small parameter. 

After introducing the column of the perturbation functions x = (^i '^)^ (upper index 
T everywhere sign the transposition) the effective action for the perturbations u{r),v{r) 
in the leading non- vanishing order ( e^) has the form 

oo 

S = j ((X+)'(X)'+P(r)(x+)'(^^2)X + X+V^(r)x-^VX) dr + u\ (13) 



where a2 is the corresponding Pauli matrix. Note, we have added the term cj^ into the 
effective action p3|) in order to get the normalization condition (see below) by a formal 
variation in respect to the A = — cu^. The elements of the matrix V{r) and the function 
p{r) are expressed in terms of the background static solution /Ar(r), $Ar(r) as follows: 



V^ii(r) 



exp{$Ar} 



2r2 



Vu{r) = V^2i(r) = ^^[rexp{<l>^/2}/ 



V2r 

V22{r) - -<S>^ + -<S>^ + — , 



n] 5 



'p(r) 



, exp{«l'^/2}/^^ 
V2r 



In order to bring the eigenvalue problem to a self-adjoined form, the extended deriva- 
tive D is introduced: 

fu'-A'v 

DX = X -^^-^X=[^,^^,^ 

Air)^j^-^^MmMhd^. (14) 
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Then the effective action ()13p can be rewritten in the form: 

S= I (Dx^Dx + X^U{r)x-uj\+x) dr + 



(15) 



The elements of the matrix U (r) are expressed as: 

exp{<l>iv} 



f/ii(r) 
Uuir) 

U22{r) 



2r2 
U2i{r) : 



r/ 2 (/at - 1)^ 
JN -I -3. 



[r exp{<l>^/2}/ 



V21 



n] ) 



(16) 



The gauge transformation 

X = exp{i A{r)a2} \E'(r) , 
where \1/ = (\l/i, ^1/2)^, allows one to bring the action to the final self-adjoint form 



(17) 



00 

S = j (($+)'(^)' + ^+f/(r)^ -c<;2^+^) dr + 



where the elements of the matrix potential U{r) are 

U{r) = exp{—iA{r)a2}U{r)exp{iA{r)a2}, 



(19) 



The behavior of the matrix elements Uij i]19\i . corresponding to the static solutions with 
= 1, 2, 3, 4 nodes, is demonstrated on Fig. El 

Let us put u'^ = —A, then we get a matrix Sturm-Liouville problem as a result of 
the unconditional extremum conditions of the functional ()18|1 on the set of the variables 
(^1,^2, A): 

-^" + f/(r)^ + A^ = 0, (20) 
on the semi-axes < r < 00 with the boundary conditions: 



^i(0) = 0, ^2(0) = 0, ^i(oo)=0, ^2(00) = 0. 
and the norm condition 



00 00 
/(^) = j ^^J+^^/dr - 1 = j {^j + ^f) dr - 1 = 0. 



(21) 



(22) 



The linear stability of the static solutions is related to the spectrum of the eigenvalue 
problem pn|) - (j22I): the negative eigenvalues u"^ = —X < correspond to the unstable 
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modes. In the Ref. ^ the Calogero phase functions method j^H] has been apphed in order 
to prove the unstable eigenmodes existence: it was shown that each static regular solution 
with nodes of the YM function had exactly unstable eigenmodes in the spectrum of 
the Sturm-Liouville problem ^ - (j22I)- 

The goal of our present paper is to get numerically the unstable eigenvalues and the 
corresponding eigenfunctions of the Sturm-Liouville problem (j^UI) - (j^^ with accuracy, 
which should be sufficient for further simulation of the static regular solution decay via 
their unstable eigenmodes. 



3 Iterative scheme for the Sturm-Liouville problem 

In what follows, the elaborated continuous analogue of Newton's method (CANM) 
P7] - [30] was applied to solve our Sturm-Liouville problem. 

The first step is to replace the semi-infinite interval r G [0, oo) by the interval r G 
[0,i?oo] for the numerical solution of the problem. Here Roo ^ is the "actual infinity". 
At the same time the problem ()2()|1 - can be rewritten in the form 

£;(^,A) = -^" + f/^ + A^ = 0, (23a) 
^^(^(0), ^'(0), A) = , Gr{^{R^), ^'{Roo), A) = , (23b) 

,(»)./*.«._i./[^;„-,*^M],._i.o, (23c) 



where Gl and Gr are boundary conditions (of mixed type, in general) at the left (r = 0) 
and at the right (r = Roo) ends of the interval. 

The boundary conditions at the left side Gl are already imposed in ((211), so we just 
rewrite it as 

Gii(^(0),^'(0),A) = ^i(0) = 0, 

G2l(v1>(0),vI/'(0),A) = %{0) = 0. (24) 

The consistency of the imposed boundary condition Gn = G2L = is confirmed by the 
local solutions behavior, obtained from the series expansions of the system ()23a|l at the 
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vicinity of the origin: 



f/ii(rUo = 263(6-4)r2 + 0(/) 

Ui2{r)r^o = U2i{r)r^o = t^b^b-3)r + Oir^), 



^ + 36 (6 - 2) + - 6^(4 + 46 - fe^) + O (r^) ; 



1 



^2(r).-.o = C'2r2 + 0(/) . 

Here Ci and C2 are free parameters. The static solutions behavior © was used to obtain 
the matrix elements Uij{r) asimptotics at r — > 0. 

The boundary conditions Gr at the right end (r = i?oo) can be defined in a similar way, 
taking into account the local solutions behavior at r ^ 00. Indeed, using the asymptotic 
behavior of the static solutions 



and the matrix elements Uij{r) 



U22{r)r- 



2 6c + d 
2 Qc + d 



sin\A^) + 0{r-^), 

-;^ + ^) sin(2 A.) + 0(r-) 

cos2(Aoo) +0(r-^). 



(25) 



where the value A^o = A{r = 00) (according to (HH)), we obtain the local solutions for 
the functions given by the following series expansion: 



C2 1 + 



COS^(y4c 



C2^^^^^^ + 0(r-2) 



2^/Xr 
sin(2y4oo 



+ 0(r- 



(26) 



VXr J 2v^r 

If we restrict our considerations within the terms of order 0(r^^) e"^'', the following 
couple {GiR, G2r) of the boundary conditions at the right end r = _Roo are given by 



G'lR(^(i?oo),^'(i?oo), A) = ^l'(i?oo) + V^^l(i?oo) = 0. 
G2i?(^(i?oo),^'(i?oo), A) = %{Roo) + V^^2(i?oo) = 0. 

The problem (j^^ is equivalent to the nonlinear functional equation 



(27a) 
(27b) 



( E{z) \ 
Gl{z) 
Gr{z) 



{2t 
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where the elements z = (A, A G M, and \&i(r) belong to the set of twice different iable 
functions on the finite interval [0,-Roo], which satisfy the boundary conditions ()23b|) . 

It is supposed that an isolated solution z* = (A*, of the equation (j28|l does exist. 
Then, according to the CANM method [2Zj - [SD], a formal evolution parameter t is 
introduced so, the following evolution equation 

^ipiz{t),t) = -ip{z{t),t), 0<t<oo, (29) 
at 

is considered with the initial condition 

It is also supposed that the given initial data zq = (Aq, ^I/q) is chosen in some vicinity 
of the exact desired solution z* in a corresponding functional space. In Ref. [^01 it was 
shown that if additionally: 

1. ip{z)is a smooth function; 

2. the operator [(y9'(z)]^^ exists and is bounded in some vicinity of the point z{0) = zo; 
then the equation has a unique local solution z(t) in the vicinity of z* , and 

lim z{t) = z* . 

t—*oo 

We assume the above two conditions are fulfilled. 

Let us introduce the functions w{r,t) and fi{t) in the following way: 

, , (9^(r,t) , , „, , , , 



dt 



_ fwi 






[r,t) 


d\{t) 




dt 





(30b) 



Then the equation ()29j) in our particular example is rewritten in terms of a couple of equa- 
tions which allow one to determine the introduced functions w{r,t) and fi{t) as follows: 



-w" + {U + XI)w = -{-"^ " + {U + fil)"^) , (31a) 

j-GLi^iR^),^'iR^),\) = -Gi(vl/(i?),vl/'(i?^),A), (31b) 

-^Gni^HiR^), ^'{Roo),\) = -Gni^iR), vl/ '{R^), A) , (31c) 
dt 



/ [^i(r,t)ti;i(r,t) + ^2(r,t)w2(r,t)] dr - - 



Roo 

2 



(^f(r,t) + '^2ir,'t))dr - 1 







0, (31d) 
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where I is 2x2 unit matrix. 

In order to solve the PDE evolution problem (jHUj) . (jHT|) numerically we have used a 
standard finite-difference technique. Let us t^,t^, . . . ,t^, . . ., t° = 0, t^^^ — t'^ = r'^ is a 
given discretization of the "time" t. We assume the step < < 1 is properly specified. 
Then we use the Euler scheme in order to obtain the next + l)-th approximation to the 
exact solutions for the eigenfunction "^^^^ = \['(r, and the eigenvalues A'^^^ = X{t^^^): 



^'=+i(r) = ^i/\r) + T^w\r), (32a) 
X^+i = X'' + T\fi'' - X''). (32b) 

It is supposed that \E''^(r) and are known at the k^^ iteration step, whereas the \E'°(r) 
and A° are some properly chosen initial data. 

At the each step in evolution k the iterative corrections w^{r) are calculated by the 
following way: 

^k^^k^^k^k^ (33) 

where the functions C,^ = (Ci , Cl )^ solutions of the problem (for sake of simplicity we 
will henceforth omit the number of iterations k): 



-C 1 + {Un + A) Ci + UuC2 = ^ 1 - f/ii^i - f/i2^2, (34a) 

-C2 + {U22 + A) C2 + UuCi = ^2 - - f/i2^i, (34b) 

Ci(0) = -vl/i(0), (2(0) = -^2(0), (34c) 

([(Roo) + VACi(^oo) = -^[(Roo) - ^ ^l(i?oo), (34d) 

C'2(i?oo) + V\C2{Roo) = -^2(^00) - ^ *2(i?oo), (34e) 

whereas the functions rj^ = (r^J^, r/f)^ are solutions of the problem: 

-Vi + iUn + X)vi + Ui2m = -"^u (35a) 

-V 2 + {U22 + A) r/2 + U12V1 = -^2, (35b) 

r?i(0)=0, r/2(0) = 0, (35c) 

r/;(i?oo) + v^r/i(i?oo) = --^^i(i?oo), (35d) 

2v A 

v'^iRoo) + ^V2{Roo) = ^^2(i?oo). (35e) 



We used the finite difference method for the discretization in respect to the spatial 
variable r. The corresponding discrete linear problem (|34j) and (|35|) has been solved by a 
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matrix sweat method. After that the value of the parameter /x'^ was calculated by means 
of the expression ()31d|) 

1 - / (^1^ + ^2') dr-2 J (^iCi + ^2(2) dr 
1^= Ji • (36) 

2 / (^iT^i + ^2^72) dr 


After the k -th approximation, \E''^ and A'^ are found, the next. A; + 1 -th one, 
and A^"^^ are calculated according to the (jHS))- The iteration process goes on until the 
fulfillment of the inequality 6^ < e, where e > is a given small value, and the discrepancy 
6^ is defined by the following relation: 

6' = MX\^')\\^ max MX',^')\. 

The best convergence of the method was achieved for the following choice jHOj of the step 
A; > 0: 



min (1, r'^-i 6''-^/6'') , 6'' < 6''~\ 

(37) 

max (rO, t''-^ 5^-^/5^) , 5^ > 6''-^ , 

for some given tq which depends on the initial approximation zq. 

In order to check and compare our results, we also have solved the following non 
self-adjoint Sturm-Liouville problem 

-X" + P(r)x' + Q(r)x = -Ax, < r < 00 , (38) 

which follows from the perturbation functional (fTSj) . Here P{r) and Q{r) are 2x2 
matrices, and P(r) is an antisymmetrical one. The boundary conditions are analogous to 
the conditions pT]) . 

We used the collocation method for the discretization of appropriate linearized bound- 
ary problems at each iteration step. The accuracy of appropriate difference scheme on 
the analytical grid with the step h is 0{h'^). The choice of the iteration step is realized 
by the Ermakov-Kalitkin formula jSO]- 

4 Numerical results 

To obtain the eigenmodes of the regular static YMD solutions the following problems 
have been solved: 
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a) the static background solutions, labelled by the number of the nodes of YM 
function, which are solutions to the boundary- value problem (jH), ©; 

b) the eigenmodes, which are solutions to the self-adjoint Sturm-Liouville problem ()20|) 
- (j221) or (equivalently) to the non self-adjoint Sturm-Liouville problem 

Both the problems a) and b) were solved using the continuous analogue of Newton's 
method. 

The semi-infinity interval [0, oo) was changed to the interval [0, i?oo] for the numerical 
reasons. The influence of the "actual infinity" Roo was studied by the establishment 
method. The results for the case = 1 on the mesh exponentially condensed to the 
origin r = 0, are shown in Tabled Here NM is the number of elements of the mesh. It is 

Tab. 1: Influence of the parameter i?oo on the eigenvalue A} 



Rco 


NM 


\\ = -u^ 


50 


101 


0.0905620823346630 


100 


201 


0.0905653977396295 


200 


401 


0.0905657469757358 


400 


801 


0.0905657826471273 


800 


1601 


0.0905657861599309 


1600 


3201 


0.0905657864933593 



evident that each doubled value of R^o leads to the establishment at least one significant 
digit after decimal point. 

The one of the specific features of the problem is the rapid increasing of Roo as the 
number (which labels the corresponding background solution) grows. Table |21 shows 
sufficient values of R^o needed for calculation of the unstable eigenmodes for the back- 
ground solutions with = 1, ...,6. 

In Table 2 the corresponding first eigenvalues {A^vjiLi are shown, which correspond 
to the background solutions with A^ = 1, 2, ..6. 

All the values of b are obtained numerically by means of the comparison of the back- 
ground solutions to the boundary- value problem (jH), (0) and the local expansions (0). In 
Ref. it is shown that the main background solution parameters tend to some limit 
values as N oo. This limit solution can be characterized by the parameter boo ~ 1.518. 
This allows us to determine the first limit eigenvalue using an extrapolation of X]^ as a 
function of A^: ^ 4.22 x 10"^ (See Table EJ. 

In Table El all the eigenvalues {A'^}^^, which correspond to the background solutions 
with A^ = 1,2,3,4, are presented. First three eigenvalues for the background A^ = 4 
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Tab. 2: First eigenvalues Ajy and appropriate parameters i?oo and b. 



N 






b 


1 


9.0566 X 10- 


-2 


2 X 10 


1.043320582 


2 


7.5382 X 10- 


-2 


2 X 10^ 


1.414072399 


3 


4.9346 X 10- 


-2 


2 X 10^ 


1.500007215 


4 


4.3455 X 10- 


-2 


1 X 10^ 


1.515017863 


5 


4.2434 X 10- 


-2 


1 X 10^ 


1.517493316 


6 


4.2266 X 10- 


-2 


3 X 10^° 


1.517897653 


oo 


^ 4.22 X 10- 


-2 




^ 1.518 



solutions are also presented. For a given background solution, labelled by the number 
A^, the fast decreasing of the value AJ^ as i increases is obvious. Hence, only the main 
(minimal) eigenmodes are presented below. 



Tab. 3: Eigenvalues {Ajylili • 



N 


A]v 






X4 


1 


9.0566 X 10- 


^2 








2 


7.5382 X 10- 


-2 


2.0742 X 10"^ 






3 


4.9346 X 10- 


^2 


1.4957 X 10-^ 


1.9622 X 10-^ 




4 


4.3455 X 10- 


-2 


5.9905 X 10-5 


1.3278 X 10-^ 


~ 10-9 



The eigenfunctions \E'|(r) and '^^i''^) for = 1,2,3,4 are shown in Fig. El The 
eigenfunctions \I'f(r) and \E'2('") for = 2,3,4 are shown in Fig. 01 The eigenfunctions 
\E''^(r) and \E'2(^) = 3,4 are shown in Fig. O 
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Fig. 3: Eigenfunctions *J and ^\ for N = 1,2,3,4. 

5 Conclusions and outlook 

We have found the eigenvalues and the eigenfunctions for the matrix Sturm-Liouville 
problem which are proper unstable modes of the regular static spherically symmetric 
solutions in the coupled system of Yang-Mills-dilaton equations. The efficiency of the 
continuous analogue of Newton's method (CANM), which was used to solve the corre- 
sponding numerical problem, is demonstrated. The CANM method, being an universal 
tool for numerical solution of both boundary value and Sturm-Liouville problems, is de- 
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Fig. 4: Eigenfunctions and for N = 2,3,4. 

scribed in detail. 

The obtained unstable modes are used as initial data in the static solutions decay 
problem in a nonlinear regime. This task is in progress now and will be reported soon [31]. 
The mentioned nonlinear decay problem arises in a wide class of mathematical physics 
problems, inspired by the recent progress in the singularity formation understanding in 
nonlinear wave equations and in massless fields collapse problem. 
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Fig. 5: Eigenfunctions *f and *2 for N = 3,4. 
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A Appendix 

To facilitate the readers, the continuous analogue of Newton's method (CANM) is pre- 
sented briefly below. 
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Let us consider the nonlinear equation 

X{x) = 0, (39) 

where the nonhnear operator y = xi^) maps the Banach space X (the element x G X) 
into the Banach space Y{y G Y). Solution to the eq. (fH^ means finding such elements 
X* G X, which are mapped by the operator x into the trivial element of space Y. 

If we suppose that Xk is known approximation to the sought solution at the kth 
iteration stage, then the increment Ax^ is computed by the formula 

Axk = ip{xk), Xk+i = Xk + Axk, k = 0,1,2, ... , where xq - given element. 

The used method of iteration determines the manner of constructing the function 
■^(x). For example, in the case of Newton's method ?/'(x) = —[x'{x)]~^x{x) where x'i^) is 
a linear operator, it is the Frechet derivative of function xi^)- Fo^^ each iteration process of 
the kind mentioned above one can build a continuous analogue introducing a continuous 
parameter t, < t < oo instead the discrete variable k, {k = 0,1,2,...). Further we 
suppose the smooth dependence x = x{t) and introduce the derivative ^a;(t) instead of 
the increment Axk- In this way we obtain the differential equation 

^x(t) = -^[(x(t)], x(0) = Xo. (40) 
at 

Thus, the solving the original eq. (j39|) is realized by the solving of the above Cauchy 
problem ()40|) on the positive half-axis < t < oo. 

There have been proved numerous theorems (see [2Zj, [2111201) related to the conver- 
gence to the isolated solution x* of the continuous analogues of various iteration methods. 

As a particular case of the continuous analogue of Newton's method (CANM) we can 
present the eq. ipUl) as 

j^x[x{t)] = -x[x{t)], x{0)=Xo. (41) 
From here we obtain the first integral 

X[x{t)] = e~\{xo). 

If the function xi^) is smooth and the operator [x'(a;)]^^ is bounded in the vicinity of 
initial approximation xq, then in the same vicinity there exists an isolated root x* of the 
eq.dnni) and limt^oo a;(t) x*. 

For example, let us consider the convergence conditions of CANM for the following 
simple non-linear boundary value problem: 

xiy) = {y" + f{x,y),y{o),y{i)} = o, xe{o,i). (42) 
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Theorem 1 Let the solution of BVP l{4^ exist and can be localized. Furthermore 

(i) the function f{x,y) is smooth in some domain D; 

(ii) the boundary value problem 

v" + f;{x,y)v = 0, viO)=v{l) = 
has only a trivial solution for each smooth function y{x) G D; 

(iii) ||?/o + f{x,yo)\\ ^ e where e > is little enough, and function yo{x), smooth in 
domain D is an initial approximation of the sought solution y*{x). 

Then the system with respect to functions y{x,t) and v{x,t) 

v'L + fyi^, y)v = - [y'L + /(a;, y)] , y't = v 

with boundary conditions 

t;(0,t) = v{l,t) = 0, 

and initial condition 

y{x,0) = yo{x) 

has on half-strip s = {{x,t) : < x < 1,0 < t < oo} a unique solution, subjected to the 
condition 

lim \\y{x,t) -y*{x)\\c2[o,i] = 0. 

t — ^oo 

The most simple method for approximated integration of the problem ()40|) is the Euler 
method. Let us build the set t^, k = 0, 1,2, ... and = tk+i — tk- Then the following 
sequence of linear problems is reached: 

x{xk)vk = -xixk), Xk+i = Xk + TkVk, fc = 0,l,2,..., 

where xq is a predetermined element. When the parameter Tjt = 1, we obtain the classical 
Newton's method. 
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